function diff = int_cost1_fun( m, param)

    diff = nan(size(m));
    f = @(x) ( (param.zeta+delta_fun( x, param))./x.^(1+param.gamma_R_1));

    for i=1:length(m)
        diff(i) = integral( f, param.m_h, m(i));
    end

end
